Anomalous magneto-oscillations in two-dimensional systems 
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The frequencies of Shubnikov-de Haas oscillations have long been used to measure the unequal 
population of spin-split two-dimensional subbands in inversion asymmetric systems. We report self- 
consistent numerical calculations and experimental results which indicate that these oscillations are 
not simply related to the zero-magnetic-field spin-subband densities. 
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Spin degeneracy of the electron states in a solid is the 
combined effect of inversion symmetry in space and time. 
Both symmetry operations change the wave vector k into 
— k, but time inversion also flips the spin so that com- 
bining both we have a two-fold degeneracy of the single 
particle energies, £+(k) = £:_(k) (Ref. [|). When the 
potential through which the carriers move is inversion 
asymmetric however, the spin-orbit interaction removes 
the spin degeneracy even in the absence of an external 
magnetic field B. This 5 = spin splitting is the sub- 
ject of considerable interest because it concerns details 
of energy band structure that are important in both fun- 
damental research and electronic device applications ( 
[p[-p^ and references therein). 

The spin splitting of the single particle energies yields 
two spin subbands with different populations N±. The 
frequencies f±'^^ of longitudinal magnetoresistance os- 
cillations in small magnetic fields perpendicular to the 
plane of the system, known as Shubnikov-de Haas (SdH) 
oscillations, have often been used to measure the 

B = spin-subband densities N± following 
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Here e is the electron charge and h is Planck's constant. 
Eq. (Q) is based on a well-known semiclassical argument 
due to Onsager which relates the cyclotron motion 
at _B > with extremal cross sections of the Fermi sur- 
face at B = 0. In this paper, we test both experimen- 
tally and theoretically the validity of this procedure. We 
obtain good agreement between experimental and cal- 
culated SdH oscillations. On the other hand the cal- 
culated B — spin splitting differs substantially from 
the predictions of Eq. (j^). We will show that this dif- 
ference reflects the inapplicability of conventional Bohr- 
Sommerfeld quantization for systems with spin-orbit in- 
teraction. 

The subject of our investigation are two-dimensional 
(2D) hole systems in modulation-doped GaAs quantum 
wells (QW's). We use GaAs because high quality sam- 
ples can be grown which allow the observation of many 



SdH oscillations, and because the band-structure param- 
eters are well known 16 1 so that accurate numerical 
calculations can be performed. The crystal structure of 
GaAs is zinc blende, which is inversion asymmetric. Fur- 
thermore, a QW structure can be made asymmetric if an 
electric field E± is applied perpendicular to the plane of 
the well. Therefore, at a given 2D hole density, the B = 
spin splitting in these systems has a fixed part due to the 
bulk inversion asymmetry (BIA) , and a tunable part due 
to the structure inversion asymmetry (SIA). 

Figure ^ highlights the main findings of our paper. It 
shows the Fourier spectra of the calculated [Fig. ^a)] 
and measured [Fig. 0(b)] SdH oscillations as well as the 
expected peak positions {h/e)N± according to the cal- 
culated spin split densities N± at B = [open circles 
in Fig. |^(a)] for a 2D system with constant hole density 
N = N+ + = 3.3 ■ 10" cm-2 but with varying Ej_. 
Even around E± — 0, when we have only BIA but no 
SIA, the open circles indicate a significant spin splitting 
AA^ = A^_|_ — N^. However, the Fourier spectra in Figs. 
[l|(a ) and (b), while in good agreement with each other 
[|l7| , ^ , deviate substantially from the zero-_B spin split- 
ting: for nearly all values of E± the splitting {h/e)AN is 
significantly larger than A/ = /_^^^ — /J''^^. In particu- 
lar, near E± = only one SdH frequency is visible in both 
the measured and calculated spectra, whereas we would 
expect to obtain two frequencies |l^. In the following 
we will show how one can understand these results. We 
will briefly describe some details of our calculations and 
experiments and then discuss the physical origin of when 
and why Eq. ([|) fails. 

Our calculations are based on the methods discussed in 
Refs. [p0|-p2t. A multiband Hamiltonian [^ containing 
the bands r§, Fg and Fy is used to calculate hole states 
in the QW. It fully takes into account the spin splitting 
due to BIA and SIA. The Poisson equation is solved self- 
consistently in order to obtain the Hartree potential. We 
obtain two spin-split branches of the energy dispersion 
£±(k||) as a function of in-plane wave vector k||. How- 
ever, we do not call these branches spin-up or spin-down 
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because the eigenstates are not spin polarized, i.e., they 
contain equal contributions of up and down spinor com- 
ponents. (This reflects the fact that for B = the system 
has a vanishing magnetic moment.) From £±(k||) we ob- 
tain the population N± of these branches [ pl| . 

For the calculation of SdH oscillations we use the very 
same Hamiltonian j23j discussed above so that the re- 
sults for B = and B > are directly comparable. We 
introduce the magnetic field by replacing the in-plane 
wave- vector components with Landau raising and lower- 
ing operators in the usual way |^,^ . From the Landau 
fan chart, using a Gaussian broadening, we obtain the 
oscillatory density of states at the Fermi energy which is 
directly related to the electrical conductivity In or- 
der to match the experimental situation the Fourier spec- 
tra in Fig. |l|(a) were calculated for B between 0.20 and 
0.85 T {B-^ between 1.17 and 5.0 T'^). We note that 
the positions of the peaks in the Fourier spectra in Fig. 
|l|(a) depend only on the Landau fan chart as determined 
by the multiband Hamiltonian . A single peak in the 
Fourier spectrum corresponds to the situation that at the 
Fermi energy the spacing between Zeeman-split Landau 
levels is a fraction a of the spacing between Landau lev- 
els with adjacent Landau quantum numbers n and n+1, 
with a constant a independent of B. 

For measurements, we use Si modulation doped GaAs 
QW's grown by molecular beam epitaxy on the (113)A 
surface of an undopcd GaAs substrate. The well width of 
the sample in Fig. |^ is 200 A. Photolithography is used to 
pattern Hall bars for resistivity measurements. The sam- 
ples have metal front and back gates that control both 
the 2D hole density and E±. Measurements are done at 
a temperature of 25 mK. In order to vary Ej_ while main- 
taining constant density we first set the front gate (Vfg) 
and back gate (Vbg) biases and measure the resistivities 
as a function of B. The total 2D hole density N is de- 
duced from the Hall coefficient. Then, at small B, Vfg is 
increased and the change in the hole density is measured. 
Vbg is then reduced to recover the original density. This 
procedure changes E± while maintaining the same den- 
sity to within 3%, and allows calculation of the change 
in E± from the way the gates affect the density. 

In Fig. [l](b) we show the Fourier spectra for the 
measured magnetoresistance oscillations. Keeping in 
mind that we may not expect a strict one-to-one cor- 
respondence between the oscillatory density of states at 
the Fermi energy [Fig. 11(a)] and the magnetoresistance 
oscillations [Fig. 0(b)] the agreement is very satisfac- 
tory. However, these experimental and theoretical re- 
sults indicate a surprising discrepancy between f^'^^ and 
{h/e)N±. In the following we will discuss possible expla- 
nations of these results. 

The common interpretation [Q of SdH oscillations in 
the presence of inversion asymmetry is based on the in- 
tuitive idea that for small B the Landau levels can be 
partitioned into two sets which can be labeled by the two 



spin subbands. Each set gives rise to an SdH frequency 
which is related to the population of the respective spin 
subband according to Eq. (|l|). However, a comparison be- 
tween the (partially) spin polarized eigenstates at i? > 
and the unpolarized eigenstates at B = shows that in 
general such a partitioning of the Landau levels is not 
possible. This reflects the fact that the orbital motion of 
up and down spinor components is coupled in the pres- 
ence of spin-orbit interaction, i.e., it cannot be analyzed 
seperately. 

For many years, anomalous magneto-oscillations have 
been explained by means of magnetic breakdown [ p6| . 
In a sufficiently strong magnetic field B electrons can 
tunnel from an orbit on one part of the Fermi surface 
to an orbit on another, separated from the first by a 
small energy gap. The tunneling probability was found 
to be proportional to exp{—Bo/B), with a breakdown 
field Bq, similar to Zener tunneling [^. This brings into 
existence new orbits which, when quantized, correspond 
to additional peaks in the Fourier spectrum of the SdH 
oscillations. However, if the anomaly of the SdH oscil- 
lations reported in Fig. |l| were due to magnetic break- 
down, for E± = we would expect several frequencies 
f S'iH with different values rather than the observed sin- 
gle frequency. In a simple, semiclassical picture a single 
frequency could be explained by two equivalent orbits in 
kjj space as sketched in Fig. ^ However, the latter would 
imply that the tunneling probabilities at the junctions ji 
and j2 are equal to one (and thus independent of B). We 
remark that de Andrada e Silva et al. [|l3| studied anoma- 
lous magneto-oscillations for spin-split electrons in a 2D 
system. Their semiclassical analysis based on magnetic 
breakdown failed to predict Bq by up to a factor of three 
and AN by up to 17% (see Table HI in Ref. 

In order to understand the deviation from Eq. (0) 
visible in Fig. |^ we need to look more closely at On- 
sager's semiclassical argument |l^ which is underlying 
Eq. (|l|). It is based on Bohr-Sommerfeld quantization of 
the semiclassical motion of Bloch electrons, which is valid 
for large quantum numbers. However, spin is an inher- 
ently quantum mechanical effect, for which the semiclas- 
sical regime of large quantum numbers is not meaningful. 
Therefore Bohr-Sommerfeld quantization cannot be car- 
ried through in the usual way for systems with spin-orbit 
interaction. In a semiclassical analysis of such systems 
we have to keep spin as a discrete degree of freedom so 
that the motion in phase space becomes a multicompo- 
nent vector field 



[2^, i.e., the motion along the spin- 
split branches of the energy surface is coupled with each 
other and cannot be analyzed seperately. In this prob- 
lem geometric phases like Berry's phase enter in an 
important way which makes the semiclassical analysis of 
the motion of a particle with spin much more intricate 
than the conventional Bohr-Sommerfeld quantization. 

One may ask whether we can combine the older idea 
of magnetic breakdown with the more recent ideas on 
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Bohr-Sommerfeld quantization in the presence of spin- 
orbit interaction. Within the semiclassical theory of Ref. 

spin-flip transitions may occur at the so-called mode- 
conversion points which are points of spin degeneracy in 
phase space. Clearly these points are related to magnetic 
breakdown. However, mode-conversion points introduce 
additional complications in the theory of Ref. so that 
this theory is not applicable in the vicinity of such points. 

Clearly we can circumvent the complications of the 
semiclassical theory by doing fully quantum mechanical 
calculations as outlined above. We have performed exten- 
sive calculations and further experiments which confirm 
that the results reported here are quite common for 2D 
systems. In Ref. js) spin splitting of holes was analyzed 
for two GaAs QW's which had only a front gate. Here 
V[g changes both the total density N — + iV_ in the 
well as well as the asymmetry of the confining potential. 
For these QW's we obtain excellent agreement between 
the measured and calculated frequencies f±'^^ versus TV 
including the observation of a single SdH frequency near 
N = 3.8 • 10^^ cm^^ when the QW becomes symmetric. 
However, there is again a significant discrepancy between 
A/ and {h/e)AN. 

Our results apply to other HI-V and II- VI semicon- 
ductors whose band structures are similar to GaAs in the 
vicinity of the fundamental gap | [T^ . Our calculations in- 
dicate that the deviations from Eq. are related to the 
anisotropic terms in the Hamiltonian. If the Hamiltonian 
is axially symmetric Eq. (|l|) is fulfilled. This is consistent 
with the semiclassical analysis of spin-orbit interaction in 
Ref. where it was found that in three dimensions no 
Berry's phase occurs for spherically symmetric problems. 
We note that for holes in 2D systems the anisotropy of 
£±(k|[) is always very pronounced It is also a well- 
known feature of the Hamiltonian for electrons, in par- 
ticular for semiconductors with a larger gap 0. Up to 
now most experiments have analyzed spin splitting and 
SdH oscillations for 2D electron systems ||-||] . To lowest 
order in k the SIA induced spin splitting in these sys- 
tems is given by the so-called Rashba term which 
has axial symmetry. For this particular case it can be 
shown analytically that Eq. (0) is fulfilled. 

For different crystallographic growth directions spin 
splitting and SdH oscillations behave rather differently. 
Moreover, these quantities depend sensitively on the to- 
tal 2D hole density iV = iV+ -t- iV_ in the weU. In Fig. | 
we have plotted the calculated SdH Fourier spectra ver- 
sus E± for a GaAs QW with growth direction [110] and 
N = 3.0- 10" cm-2 [Fig. |(a)] and TV = 3.3 • 10" cm^^ 
[Fig. ||(b)]. Open circles mark the expected peak po- 
sitions {h/e)N± according to the spin splitting N± at 
B = [|l^,|l^. Again, the peak positions in the Fourier 
spectra differ considerably from the expected positions 
{h/e)N±. Close to E± = there is only one peak 
at {h/2e)N. Around E± = 1.0 kV/cm we have two 
peaks, but at even larger fields E± the central peak at 



{h/2e)N shows up again. At E± « 2.25 kV/cm we have a 
triple peak structure consisting of a broad central peak at 
{h/2e)N and two side peaks at approximately {h/e)N±. 
In Fig. we have a significantly smaller linewidth than 
in Fig. nf. Basically, this is due to the fact that for the 
Fourier transforms shown in Fig. ^ we used a significantly 
larger interval of (10.0 T~^ as compared with 3.83 
T~^) in order to resolve the much smaller splitting for 
growth direction [110]. We note that for E± = the SdH 
oscillations are perfectly regular over this large range of 
B^^ with just one frequency, which makes it rather un- 
likely that the discrepancies between A/ and {h/e)AN 
could be caused by a i? dependent rearrangement of holes 
between the Landau levels. 

Similar results like those shown in Figs. ^ and ^ have 
been obtained also for growth direction [001], but the 
spectra were more complicated with, e.g., several SdH 
frequencies for E± = 0. Our calculations for holes are 
based on the fairly complex multiband Hamiltonian of 
Ref. [Q. We obtained qualitatively the same results by 
analyzing the simpler 2x2 Hamiltonian of Ref. [|l^ . How- 
ever, this model is appropriate for electrons in large-gap 
semiconductors, where spin splitting is rather small, so 
that it is more difhcult to observe these effects experi- 
mentally. 

In summary, we have both measured and calculated 
the SdH oscillations of 2D hole systems in GaAs QW's. 
As opposed to the predictions of a semiclassical argument 
due to Onsager, we conclude that the B = spin split- 
ting is not simply related to the SdH oscillations at low 
magnetic fields [Q. This is explained by the inapplica- 
bility of conventional Bohr-Sommerfeld quantization for 
systems with spin-orbit interaction. 
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FIG. 1. Fourier spectra of the (a) calculated and 
(b) measured SdH oscillations versus magnetic field B for 
different values of electric field E±_ for a 200 A wide 
GaAs-Alo.3Gao.7As QW with growth direction [113] and 2D 
hole density A*' = 3.3 ■ 10^^ cm~^. The open circles show the 
expected Fourier transform peak positions {h/e)N± according 
to the calculated spin splitting N± at B = 0. 



FIG. 2. Qualitative sketch of the spin-split Fermi contours 
in k|| space for a QW with growth direction [113] (solid lines). 
In a simple semiclassical picture the observation of a single 
peak near E± = in the Fourier spectra of Fig. ^ can be 
explained by trajectories in ky space which follow the dashed 
lines at the junctions ji and j2. 
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FIG. 3. Calculated Fourier spectra of the SdH oscilla- 
tions versus magnetic field B for different values of electric 
field E± for a 150 A wide GaAs-Alo.sGao.s As QW with crys- 
tallographic growth direction [110] and 2D hole densities (a) 
N = 3.0 ■ 10" cm-2 and (b) = 3.3 ■ 10" cm'^. The 
open circles show the expected Fourier transform peak posi- 
tions {h/e)N± according to the calculated spin splitting N± 
at B = 0. 



